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ABSTRACT 

We present MIPS observations of the cluster A3266. About 100 spectroscopic cluster members have 
been detected at 24 /*m. The IR luminosity function (LF) in A3266 is very similar to that in the Coma 
cluster down to the detection limit L[r ~ 10 43 ergs s _1 , suggesting a universal form of the bright end 
IR LF for local rich clusters with M ~ 1O 15 M . The shape of the bright end of the A3266-Coma 
composite IR LF is not significantly different from that of nearby field galaxies, but the fraction of 
IR-bright galaxies (SFR > 0.2M Q yr _1 ) in both clusters increases with cluster-centric radius. The 
decrease of the blue galaxy fraction toward the high density cores only accounts for part of the trend; 
the fraction of red galaxies with moderate SFRs (0.2 M Q yr" 1 < SFR < 1 M Q yr _1 ) also decreases with 
increasing galaxy density. These results suggest that for the IR bright galaxies, nearby rich clusters 
are distinguished from the field by a lower star-forming galaxy fraction, but not by a change in L* IR . 
The composite IR LF of Coma and A3266 shows strong evolution when compared with the composite 
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IR LF of two z ~ 0.8 clusters, MS 1054 and RX J0152, with L* IR oc (1 + zf 

. This L* IR evolution is indistinguishable from that in the field, and the evolution is stronger, 
but still consistent with that in the field. The similarity of the evolution of bright-end IR LF in very 
different cluster and field environments suggests either this evolution is driven by the mechanism that 
works in both environments, or clusters continually replenish their star-forming galaxies from the field, 
yielding an evolution in the IR LF that is similar to the field. The mass-normalized integrated star 
formation rates (SFRs) of clusters within 0.5i?20o also evolve strongly with redshift, as (1 + z) 5 3 . 
Subject headings: galaxies: clusters: individual (Abell 3266) — galaxies: luminosity function — in- 
frared: galaxies 
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1. INTRODUCTION 

Galaxy clusters are unique laboratories to study the 
environmental effects on galaxy evolution. A high den- 
sity environment could potentially alter the star for- 
mation (SF) in galaxies through many different phys- 
ical mechanisms, e.g., ram pressure stripping, galaxy- 
galaxy interactions, a nd strangulation of the gas reser- 
voir of galaxies (see iBoselli fc Gavazzil I2006L and ref- 
erences therein). As an extreme case of high den- 
sity environment, clusters show clear evidence of less 
star formation co mpared to the fie l d at both local and 
higher redshifts (iKennicuttJ Il983t iBalogh et alJ 119971 : 
lHashimoto et al.l Qj)98; Po ggianti et alJ 11999ft . suggest- 
ing environmental mechanisms p lay a role in transform- 
ing th eir star forming properties (|Christlein &; Zabludofj 
2005). However, it is still not clear which mecha- 
nism, or mechanisms, operate and where the trans- 
formation happens in the first place. Is it an end 
produ ct of galaxy 'preprocessing' in the group environ- 
ment dZabludoff et al.lfl9 96: Zabl udoff k Mulchaevlfl998l : 



iLewis et al.l 120021 : iGomez et al.ll2003| ) or is it caused by 
some mechanism specifically related to the cluster en- 
vironment (c.f., iBerrier et al.ll20(M) ? Quantitative and 
unbiased star formation measurements in clusters and 
their comparison to field galaxies are keys to answering 
these questions. 

The luminosity function (LF) provides such a quan- 
titative tool to study galaxy properties. It has been 
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widely used to compare galaxies in clusters with those 
in the field at many different wavelengths. How- 
ever, some of the results are controversial. For ex- 
ample, some studies claim that cluster LFs show lit- 
tle variation across a wide ran ge of cluster prop erties 
(ICollesslll989HRauzv et al.lll998tlDe Propris et atll2003l : 
iPopesso et al.ll2006l ). while others find t hem to depend on 
clust er richness, Bautz-Morgan type dBautz &: Morgan! 
1970fl, or distance from the cluster center (Ipresslerlll978 



Garill i~et al.lll999t lLopez-Cruz et al.lll997l : lHansen et al.l 
20051 iBarkhouse et al.l l2007h . Another question is 

whether cluster LFs diff er from those in the fi eld. De 

Propris et al . (199 8), IChristlein fc Zabludofll (|2003h . 
ICortese et all (|2003ft and iBai et al.l (|2006h found clus- 
ter LFs to be indistinguishable from field LFs, but oth- 
ers suggest that they have brighter chara cteristic mag- 
nitudes and different faint end slopes dValotto et al 



| 1997tlGoto et al.ll2002trYagi et al.ll2002T:lDe Propris et al 
120031) . 

One of the major reasons for these apparent contradic- 
tions is that these studies are based on different wave- 
length ranges and therefore probe different properties of 
the cluster galaxy population. The optical and near- 
IR luminosities of galaxies measure the old stellar mass, 
while the blue and UV band luminosities are more sen- 
sitive to the young stellar population. Since the LFs of 
differ ent types of galaxies are different (jMadgwic k et al.1 
2002), the LFs based on different wavelengths are bi- 
ased toward different s ub-p opulations of cluster galax- 
ies. iDe Propris et al.l (|2003l ) found that although the 
B-band LF of early-type galaxies in clusters is differ- 
ent from that of the field, the B-band LF of late-type 
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galaxies in clusters is very similar to that of the field. 
A similar result was also found in the R-band for star 
formin g and quiescent g alaxie s bv lChristlein fc Zabludoflj 
IpOOlh . iCortese et all (120051) found that the UV LF of 
star-forming galaxies in a nearby cluster is consistent 
with the field and that the difference between the cluster 
and field LFs is due to the quiescent galaxies. 

Although UV and B-band luminosities are sensitive to 
star-forming galaxies, they can be heavily attenuated by 
dust and therefore bias the LF. In contrast, the total IR 
luminosity (Lir, 8-1000 /im) shows a good correlation 
with the star f ormation rate (S F R) of galaxies un affected 
by extinction (|Kennicuttlll998l ). IBai et alj (|2006f ) studied 
the IR LF of the Coma cluster and concluded that its 
shape shows no significant difference from that of the 
field IR LF. However, this result is drawn from the study 
of one rich cluster and its validity in other clusters needs 
to be tested. 

The study of the evolution of LFs with redshift can 
provide crucial constraints on models for the formation 
and evolution of galaxies. In particular, the evolution 
of the IR LFs will help us understand the star formation 
history. By comparing Coma with two high redshift clus- 
t ers, MS 1054-03 and RX J0152 (jMarcillac et al.l 120071) 
, IBai et al.l (|2007l ) found a strong evolution of the IR 
LFs for luminous galaxies in both density and luminosity 
(Lj R ) from z ~ to ~ 0.8. The evolution rate is con- 
sistent with that found in field IR LFs and suggests that 
this evolution pattern is probably not caused by clus- 
ter environmental mechanisms, but by processes similar 
to those that regulate the field IR LF evolution. How- 
ever, because the comparison was based on only one low 
redshift cluster, Coma, there are possible systematic un- 
certainties in this evolution trend, especially if the IR 
LFs of local rich clusters have a large cluster-to-cluster 
variation. It has become clear that the Coma Cluster, 
once thought to be the archetype of a well-relaxed clus- 
ter, has many small substructures dMellier et alj 119881 : 
IWhite et al] 119931: iNeumann et~al] I2003D . Although we 
did not find any significant difference between the IR LF 
of the infalling NGC 4839 group and that elsewhere in 
the Coma cluster, it is still not clear how these substruc- 
tures may affect the overall IR LF. To study how different 
dynamical states may affect the IR LF and to quantify 
cluster-to-cluster variations, we need to study more low 
redshift clusters. 

To test the universality of the IR LF in local clus- 
ters and provide a reliable baseline for both the cluster- 
field comparison and for evolution in clusters, we present 
the IR LF of A3266 in this paper. A3266 is a nearby 
rich cluster (z = 0.06) with a mass similar to Coma 
(~ 10 15 M Q ). It has a large sample of spectroscopi- 
cally confirmed cluster members (more than 300), mak- 
ing it possible to obtain a IR LF complete down to 
Lir ~ 10 43 ergs s _1 . It is very bright in the X- 
ray [Lx ~ 10 45 ergs s~ 4 ). Its X-ray morphology and 
temperature map, as well as a comparison of the core 
and overall velocity dispersions, all suggest a recent ma- 
jor merger dQuintana et alj|1996t iHenriksen et al.l F2000; 
ISauvageot et al.l 120051 : iFinoguenov et al.ll2006l ). In §2, 
we present the data analysis for this cluster. In §3, we 
calculate the IR luminosities and SFRs for the cluster 
members. In §4, we determine the IR LF for A3266, 
compare it with Coma's, and look at the dependence 



of the IR galaxy properties with cluster-centric radius. 
In §5, we discuss the implications of the results and 
study the evolution of the integrated SFRs. In §6, we 
summarize the main conclusions. Throughout this pa- 
per, we assume a ACDM cosmology with parameter set 
(fc,n 0> Ao) = (0.7,0.3,0.7). 

2. OBSERVATION AND DATA REDUCTION 
2.1. MIPS Observations and Source Extraction 

Th e A3266 cluster was observed by MIPS (|Rieke et all 
2004) on Spitzer in medium scan map mode (GTO pro- 
gram # 83, PI G. Rieke) on Jun 28, 2005 at 24, 70, 160 
/xm simultaneously. A rectangular 45' x 60' region was 
mapped. The data were processed with the MIPS Data 
Analysis Tool (DAT version 3.02; iGordon et all l2005h 
and array-averaged background subtraction was applied 
to improve the images. The final mosaics have exposure 
times of ~ 80, ~ 40 and ~ 8 s pixel -1 at 24, 70 and 160 
/im, respectively. The cluster resides in a region with 
low IR background (~ IhMJy/sr). The online SSC tool 
SENS_PET estimates that the 3 a sensitivity of these 
observations are 0.25, 15 and 150 mJy at the 24, 70 and 
160 /im bands respectively. 

We used SExtractor to extract sources and measure 
their photometric values. Before source detection, the 
images were filtered with Gaussian functions with the full 
width at half maximum (FWHM) matching the FWHM 
of the point spread function of the images to help improve 
the detection of faint sources. The flux of each source 
(flux_auto) was computed in an adaptive Kron aperture. 

2.2. Spectroscopic Data 

To derive a spectroscopically-confirmed cluster mem- 
ber list, we retrieved galaxies with 0.047 < z < 0.072 
in the cluster field from the NASA/IPAC Extragalactic 
Database (NED). This redshift range corresponds to ±3 
times the line-of-sight velocit y dispersion of the cluster 
(|Christlein fe Zabludoffl l2003f) . We retrieved 298 clus- 
ter members in the central 45' x 60' region observed by 
MIPS. 

The simple method we used here for member selection 
is consistent with more sophisticated methods that com- 
bine radial velocity and position information. In Fig. [TJ 
we plot the relative line-of-sight velocity (V — (V), (V) 
is the cluster mean) versus projected distance from the 
cluster center for all the members we select. We also plot 
the maximum line-of-sight veloci ty as a function of the 
projec ted distance as proposed bv lden Hartog fc Katgertl 
(|1996f h All the members fall within, or very close to, 
this stringent criterion of member selection. We con- 
clude that the simple velocity selection in the cluster 
region studied here works efficiently in excluding fore- 
ground/background galaxies. 

The spectroscopic data in NED come from several dif- 
ferent surveys and are not homogeneous. Therefore, we 
made us e of a J?-band photometric c atalog in the clus- 
ter field (jChristlein fe Z abludofi 2003) to define the com- 
pleteness of these data. We obtained all the galaxies with 
rcdshifts in the same region from NED and compared 
their number with the total number of galaxies in the 
photometric catalog. In Fig. [5J we show these ratios as a 
function of R magnitude. As can be seen in the plot, the 
spectroscopic data have a high degree of completeness at 
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R < 17 mag, but the completeness drops rapidly to less 
than 50% at R > 17.5 mag. 

This completeness function does not change much 
across the whole survey region, although the data are 
slightly more complete in the cluster center. When we 
calculate the IR LF of the cluster members in the fol- 
lowing section, we will use the inverse of this function as 
weighting factors to correct for the incompleteness of the 
spectroscopic survey. 

3. ANALYSIS 

3.1. Correlating the IR emission with Cluster Members 

We correlated the spectroscopically confirmed cluster 
members with the nearest 24, 70 and 160 /jm sources 
projected within 6", 10" and 15", respectively. These 
matching radii take into consideration the large FWHM 
of the images, 6", 18" and 40" at 24, 70 and 160 /mi, 
respectively. At 24 /mi, the rather large matching ra- 
dius relative to the small astrometric uncertainties (< 1") 
also accounts for the possible physical displacement be- 
tween the optical and 24 /im brightness centroids. None 
of the 24 /im sources above the 3a detection limit has 
multiple matches with the spectroscopic cluster mem- 
bers. Even when matching 24 /im sources with all the 
extended sources in the photometric catalog, only 3% of 
the mid-IR sources have more than one optical counter- 
part. This demonstrates that blending is not an issue in 
correlating the IR emission with these cluster members. 
Altogether, we found 109, 48 and 15 cluster members 
with 24, 70 and 160 /im emission, respectively, above the 
3cr detection limits. 

3.2. Deducing Total IR Luminosity from 2J^ [im 
Emission 

The total IR luminosity (A = 8 — 1000/im) of a star- 
forming galaxy is proportional to its SFR and can be 
used as a SFR indicator (e.g., iKennicu ttl UflM) . How- 
ever, because direct measurement of the total IR lumi- 
nosity is only possible for a limited number of galax- 
ies, many SFR studies base their calculation on the to- 
tal IR luminosity extrapolated from a single band IR 
measurement. The mid-IR continuum emission, e.g., in 
the Infrared Astronomical Satellite (IRAS) 12 /im and 
Spitzer 24 /mi bands, has been sh own to correlate with 
the total IR luminosity very well dTakeuchi et alj 120051 : 
lAlonso-Herrero et al.|[200l ICalzetti et al.ll2007f h 

For the sake of continuity with our previous work, 
we base our estimates of the total IR luminosities on 
a sample of s tar- fo rming galaxy SEDs developed by 
iDale fc Heloul ()2002l ). These SEDs are based on IRAS 
and Infrared Space Observatory (ISO) observations of 69 
normal star-forming galaxies and have been well cali- 
brated in the far-infrared and submillimcter bands. They 
are luminosity-dependent and their total IR luminosities 
are monotonically correlated with the 24 /mi luminosity. 
We find the IR luminosities of the cluster members from 
their 24 /im flux by interpolating the Lib — fi^ relation 
given by the template SEDs (see Table Q]). The deduced 
total IR luminosities of the cluster members above the 
3<t 24 /im detection level range from 4 x 10 42 to 3.6 x 10 44 
ergs s . 

To test if the total IR luminosities deduced from the 
single 24 /tm band are consistent with the emission at 



longer wavelengths, we compare the 70 vs. 24 /im flux 
ratio (/70//24) of the SED templates to the flux ratio of 
cluster galaxies detected at both wavelengths. Because 
of the large uncertainties in the 70 /im flux measurement 
of faint sources, we only include galaxies detected at 70 
/im with a significance > 6cr. For 43 galaxies with secure 
24 and 70 /im flux measurements, the average /70//24 
is 13, only slightly larger than the /V0//24 ~ 12 given 
by the SEDs. We also determine the total IR luminosi- 
ties from the observed 70 /im flux by interpolating the 
Lib. ~ /70 relation given by the template SEDs. The total 
IR luminosities deduced in this way are consistent with 
those deduced from the 24 /im flux, with a small disper- 
sion of ~ 0.1 dex. This consistency confirms the validity 
of using the 24 /tm flux alone for the total IR luminosity 
determination. Therefore, the following studies are all 
based on the 24 /im data only. 

3.3. Comparison of SFRs Deduced from Different 
Methods 

Once we have the total IR luminosities of the 
galaxies, we calculate SFRs using the standard for- 
mula SFRlMyy r- 1 ] = 4.5 x lQ- 44 L f R \erqs s' 1 ] from 
Kennicuttl (119981). Recent studies (lAlonso-Herrero et al.l 
2006t iPerez-Gonzalez et~aTI l2006t ICalzetti et al.1 l2007t 
Rieke et al. 2008) have shown that the 24 /im emis- 
sion correlates tightly with SFR and suggest using the 
24 /im emission directly as a SFR indicator to mini- 
mize the uncertain ties associat e d wit h deriving the to- 
tal IR luminosity. iRieke et all (|2008t) find that the 24 
/im flux provides accurate SFR estimates over a large in- 
frared luminosity range. For normal star-forming galax- 
ies with 1O 1O L0 < Lib < lO 11 -^©, they provide a rela- 
tion SFR(M Q i/r~ :L ) = 7.8 x 10- 10 [L 2 4/im/£o] accurate 
within 0.2 dex. We calculate the SFRs with this formula 
and compare them with the SFRs that we derived using 
our SED fitting routine. We find good consistency be- 
tween these two SFR estimates, with a systematic offset 
of ~ 0.3 dex. The systematic offset arises primarily from 
th e different stellar initi al mass functions (IMFs) used 
in iKennicuttl (j 19981 ) and IRieke et~aTl (j2008f l. The offset 
is reduced to ~ 0.1 dex once allowance is made for the 
difference introduced by the IMFs. We will use the SFR 
obtained from the Dale & Helou templates in this paper. 
This approach facilitates comparison with previous work 
since systematic errors in different SFR derivations are 
avoided. 

For the luminosity ranges of interest in this work, the 
UV-emitting (unobscured) component of star formation 
is a small frac tion of the infrared-emitting one (e.g., 
iBai et all 12001 . so we have not tried to correct for it. 
There is no reason to expect such a correction to evolve 
strongly with redshift over the z < 0.8 range relevant to 
this paper. The infrared measurements as used here pro- 
vide a single, uniform, and accurate metric for the SFR 
over the entire redshift and luminosity range of interest. 

4. RESULTS 

4.1. IR Luminosity Function of A3266 

We calculate the IR LF over the 16 Mpc 2 surveyed 
area in A3266. The number of IR cluster members in 
each luminosity bin is divided by the surveyed area to 
obtain the surface density of IR cluster galaxies in that 
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bin. The results are shown as the open squares in Fig. [3j 
the vertical clotted line is the 24 (im detection limit. To 
account for the incompleteness of the spectroscopic data, 
we use the inverse of the completeness function (shown 
in Fig. [2]) as weighting factors when we calculate the 
number counts in each luminosity bin. 

The filled squares are the results of the incompleteness- 
corrected IR LF. The correction raises the LF slightly at 
the faint end, but there is still a drop of galaxy num- 
ber density at the lowest luminosity. The drop is caused 
by the limits in the sensitivity of the spectroscopic sur- 
vey. If we calculate the median fn/f 24 ratio of the cluster 
members with 17 < R < 17.75, where the spectroscopy is 
around 50% complete, and assume that the fainter galax- 
ies have similar /r//24 ratios, the median Lj R is about 
10 42 - 7 ergs s" 1 for galaxies at R = 19. Although the 
IR data are complete to this luminosity level, only about 
5% of the cluster members have spectroscopy (Fig. [2]). 
Thus, we start to lose IR galaxy number counts due to 
the detection limit of the spectroscopic survey, explaining 
the drop in the lowest luminosity bin. Another problem 
of the simple strategy we adopt to correct the incom- 
pleteness is the color bias between the cluster and field 
galaxies. We calculate the spectroscopic completeness in 
the R band, but at least 70% o f the spectroscopic tar- 
gets are selected in the B band (jChristlein fc Zabludofj 
l2003| h Because field galaxies are generally bluer than 
cluster galaxies, this color bias will cause an overesti- 
mate of the spectroscopic completeness in the R band 
for cluster galaxies and therefore an underestimate of the 
number density of IR cluster galaxies. We expect that 
this color bias does not affect the data points brighter 
than Ljjj w 10 43 ergs s" 1 , because most of these IR 
bright galaxies are blue in color too. However, many of 
the IR galaxies with Li R < 10 43 ergs s _1 are in the 
red sequence (as suggested by the Coma cluster data), 
and the spectroscopic incompleteness correction that we 
make may not be adequate there. 

Along with the I R LFs of A3266 , we plot the IR LF of 
the Coma cluster (|Bai et all 1200(1 in Fig. H The Coma 
IR LF is obtained with a similar data set as in A3266. 
However, iBai et all (|2006f ) use slightly different SEDs to 
deduce the total IR luminosity from the 24 /jm flux den- 
sity. We update the total IR luminosities of the Coma 
galaxies using the same SEDs that we use in this paper. 
The resulting Coma IR LF is only slightly changed. 

From the plot, we can see that the IR LF of A3266 
is very similar to that of the Coma cluster in the lumi- 
nosity region above the completeness limit. For further 
comparison, we fit the incompleteness-corrected IR LF of 
A3266 with a Schechter function (4> cx (L/ 'L*y a e- L I L \ 
ISchechterl I1976T) . We adop t the chi-sq uare minimization 
method used in IBai et all (|2006l . l2007h to find the best 
fitted parameters. The method incorporates the non- 
detection of galaxies beyond the brightest bin into the 
chi-square calculation. Because we do not have many 
data points to constrain the faint end slope, we have 
fixed it at a = 1.41, the value given by the Coma IR 
LF. We exclude the faintest data point, which is affected 
by the inadequate incompleteness correction, from the 
fitting. The resulting best fitted parameters are: 

a = 1.41 (fixed), log(L* IR /L Q ) = 10.49±°;fr (1) 

The L*j „ value is the same as the value found in the Coma 



cluster (\og{L* IR / Lq) = 10.49 + g;^. lBai et alEioM) . The 
close match in both LF shape and L* IR demonstrates that 
the bright-end IR LFs are virtually identical. A similar 
result has been obt ained previously for two high redshift 
clusters (z = 0.83: IBai et al.ll2007l ). suggesting that the 
variation of the LF from one massive cluster to another 
at the same redshift is likely to be small. 

4.2. Composite IR Luminosity Function 

The similarity of the IR LFs of the Coma cluster and 
A3266 allows us to obtain a composite LF for these two 
clusters. For this purpose, we include only the IR clus- 
ter members with Lj R > 10 42 ' 7 ergs s _1_1 . Above this 
limit, the IR surveys in these two clusters are both nearly 
complete, and the spectroscopic incompleteness correc- 
tion in A3266 works well. The composite IR LF is shown 
in Fig. [H The best fitted parameters of the Schechter 
function are: 

a = 1.41 (fixed), log(L* IR /L Q ) = 10.50±^ 2 . (2) 

The L* IR of the composite LF is well within the la error 
of the values obtained for the Coma and A3266 clusters 

i ndividually . 

IBai et al l (|2007ft found the IR LF of a high redshift 
cluster, MS 1054-03 (z = 0.83), shows strong evolution 
compared with the IR LF of the Coma cluster. The best 
fitting Schechter function of MS 1054-03's IR LF has a 
Lj R about an order of magnitude larger than that found 
in the Coma cluster. To further test this trend, here we 
combine the data from MS 1054-03 wit h the data from 
anoth er high redshift cluster, RX J0152 (Marcillac et al.l 
120071) . and obtain a composite IR LF. The total IR lu- 
minosities are deduced in the same way to minimize the 
difference caused by systematic errors. The composite 
IR LF of the high-z clusters has logL* IR /L G = 11. 35±£;i 4 
with a fixed a = 1.41, very close to the value obtained 
for the IR LF of MS 1054-03 alone. The normalization 
is slightly smaller than the IR LF of MS 1054-03, but 
is still much higher than the composite LF of the low- 
z clusters. The L* IR value increases by almost an or- 
der of magnitude from z ~ to 0.8. This result con- 
firms the strong evolution of IR LFs that we found by 
comparing only the Coma IR LF to that of MS 1054- 
03. A strong evolution is also found in the field IR LF 
from z = up to z ~ 1. iLe Floc'h et alj (|2005f) char- 
acterize the evolujjipn of IR LF in th^ Q CDF-S field as 
L* T „ oc (1 + z) 3 - 2 -°'-2 $* r „ cx (1 + z)°- 7 -° up to z ~ 1; 
iPerez-Gonzalez et al.l (|2005l ) found a similar result using 
a larger data set. In comparison, the evolution found be- 
tween the low and high-z cluster LFs can be described by 

L* IR cx (1 + z) 3 ' 2 -",^ oc (1 + z) 1 - 7 -! °. The evolution 
of the L* IR in clusters agrees with the evolution in the 
field very well. The evolution of the normalization (&*j R ) 
in clusters is stronger than that in the field, but it is con- 
sistent with the field given the error of the measurement. 
In Fig 21 we compare the evolution of IR LFs in clusters 
directly to that in the field by evolving the composite IR 
LF of Coma and A3266 to z = 0.83 using the evolution 
trend given by the field IR LF. As shown in the plot, the 
evolved LF is similar in shape to that observed in the two 
high-z clusters, b ut with a smaller n ormalization. Using 
Ha emission lines. iFTrm" et al.l (|2008f ) studied the star for- 
mation properties of a large sample of local clusters and 
compared it with several z ~ 0.75 clusters. Consistent 
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with our results, they found that the total Ha luminos- 
ity of the high-z clusters is lOx that of the local clusters 
and this evolution is comparable to the Ha luminosity 
evolution found in the field. 

It is worth noting that the difference of IR LFs between 
our high-z and low-z cluster samples is primarily due to 
evolution, not to the difference in other cluster prop- 
erties, e.g., mass. Although these four clusters are all 
among the richest clusters at their epochs, the two low-z 
clusters are more massive than the two high-z clusters. 
However, the mass of the Coma is only slightly greater 
than that of MS 1054-03 (1.4 x 10 15 M Q vs. 1.1 x 10 15 
M©). A3266 is twice as massive as Coma, but its IR LF 
is very similar to Coma's; RX J0152 is half of the mass 
of MS 1054, but ag ain, their IR LFs are not very differ- 
ent (|Bai et al.ll2007f ) . The similarity suggests that at the 
same epoch, the masses of these rich clusters do not affect 
the properties of their IR LFs substantially. Although 
some evidence of a correlation b etween star formatio n 
and cluster mass has been found (jPoggianti et al.| [2006) . 
star formation is almost inde pendent of cluste r mass for 
the mass range probed here (|Finn et al.ll2008l) . 

Another issue that could bias the comparison is the 
different aperture of the clusters within which we calcu- 
late their IR LFs. Since the cluster IR LF is the surface 
density averaged over the survey area, the different data 
cove rage in each cluster could cause differences in their 
LFs ()Popesso et al.ll2~006D . For the three rich clusters we 
study here, Coma, A3266 and MS 1054-03, the IR LFs 
are all calculated in a region about 70% of the i?200 area 
(R200 is the radius within which the mean cluster density 
is 200 times the critical density of the universe at that 
redshift). The data on RX J0152 covered the whole i? 2 oo 
area, however, i ts IR LF does n ot differ much from that 
of MS 1054-03 (|Bai et alj 120071) . In addition, in the fol- 
lowing section ( N4.3p . we show that the IR LFs over the 
majority of the volumes of the two local clusters are con- 
sistent with their overall IR LFs. This result indicates 
that the aperture difference does not affect the evolution 
we find between high-z and low-z cluster IR LFs. 

4.3. Composite IR LFs at Different Radii 

iBai et all (|2006l ) studied the IR LFs in different re- 
gions of the Coma cluster and found a flatter faint end 
slope and lack of very bright galaxies in the core re- 
gion. Although the data in A3266 are not deep enough 
to constrain the faint end slope of its LF, the composite 
LFs combining similar radial regions of these two clusters 
have better statistics at the bright end and can constrain 
variations of the bright-end LF at different cluster radii. 
For this purpose, we combine the same regions of the 
two clusters relative to i?200- -R200 is better than a fixed 
physical scale because it takes int o account the differ- 
ences in the masses of the clusters ()Popesso et al.|[2006t 
Barkhous e et al.ll2007f ). We define three regions in each 
cluster: the core region (r < 0.13i?2oo), the intermediate 
region (0.13i?2oo < r < 0.3i?2oo), and the outer region 
(0.3i?2oo < r < 0.55i?2oo)- The core region is chosen 
to m atch that define d earlier in the Coma cluster (~ 0.3 
Mpc, IBai et af. 2006) . The largest radius in the outer re- 
gion is limited by the data coverage. The three regions in 
A3266 and Coma are shown in Fig. Altogether, there 
are 8, 23 and 44 IR galaxies with Lm > 10 43 ergs s _1 
found in the core, intermediate, and outer regions, re- 



spectively. 

The resulting composite IR LFs in the three different 
regions are shown in Fig. [SJ In the core region, despite 
the high overdensity of cluster members, there are only 
8 of them with Lm > 10 ergs s _1 , and no galaxies in 
the brightest bin of the LF. Because the data points are 
so few in this region, we can not derive a reasonable L* IR 
value. Therefore, we fix both L* IR and a to the values 
of the total composite IR LF from H4.2I and only allow 
the normalization to change. The best-fitted Schechter 
function has a slightly larger normalization than that 
of the total composite IR LF. To further test whether 
those IR bright galaxies really reside in the cluster cores 
or are galaxies at larger radii being projected on to these 
regions, we compare their absolute relative line-of-sight 
velocity (|V — (V)\) with the core galaxies that are not 
IR bright (see Fig. [T|). In the A3266 core, the median 
value of \V — {V)\ for the IR bright galaxies is ~ 1600 
km s _1 , much larger than the median of non-IR bright 
galaxies (~ 800 km s _1 ). In Coma, there are only two 
galaxies with Lm > 10 43 ergs s" 1 in the core region, 
and the bright one has an absolute relative velocity about 
2.5 times larger than the median value of the rest of the 
core members. These results suggest that the few core 
galaxies with Lm > 10 43 ergs s _1 are likely either to be 
projected on the core, or to be on orbits systematically 
different from those of the other core galaxies. In cither 
case, these galaxies are probably not among the well- 
relaxed core galaxy population. 

In the intermediate and outer regions, we fit the com- 
posite LFs with a Schechter function assuming a fixed 
faint end slope as before. The difference between the in- 
termediate region (logij^ = 44.15io.24) and the outer 
region (logL* IR — 44.00^" '^ 3 ) is within the la error, and 
both L* IR values are consistent with the value for the total 
composite IR LF. Limited by the small number counts, 
we are unable to perform a statistically significant com- 
parison between the core IR LF and the LFs in other 
regions. Overall, except for the quite uncertain core IR 
LF, the IR LFs in different regions of the clusters are 
generally consistent with each other and with the overall 
composite IR LF. 

4.4. Fractions of IR-bright Galaxies 

Although the IR LFs of Coma and A3266 are al- 
most identical at the bright end, the fractions of IR- 
bright galaxies in these two clusters are slightly dif- 
ferent. We calculate the fractions for all the galaxies 
with Mr < —20.15. This magnitude corresponds to 
R = 17.0 and 14.9 in A3266 and Coma, respectively, 
both above the 60% spectroscopic completeness limits 
in their fields. In A3266, we find that 40± 3 % of cluster 
members brighter than Mr = —20.15 have Lm > 10 42 7 
ergs s^ 1 (SFR> 0.2 M Q yr _1 ) and in Coma this frac- 
tion is 31±f%. If we include the galaxies one magni- 
tude fainter (Mr < —19.15), the fractions are lower, 
35±l% and 23±l% for A3266 and Coma, respectively. 
This fainter magnitude cut includes all (except a cou- 
ple) of the galaxies with Lm, > 10 42 7 ergs s _1 in both 
clusters and it co rresponds to a galax y stellar mass of 
3 - 7 x 10 9 M Q ([Bell fc de Jongll200l . In A3266 this 
magnitude extends to the range where the spectroscopic 
data are rather incomplete and the calculation might not 
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be reliable. We note that although A3266 has a higher 
fraction, and also a higher number, of IR-bright galaxies 
than Coma, they are distributed over a larger area and 
therefore the normalization of its IR LF is still similar to 
that in Coma. 

From the velocity vs. projected radial distance plot 
(Fig. Q} and the sky map of A3266 (Fig. O, we can see 
that these IR-bright galaxies scatter around the whole 
region, but appear to be more frequent at large radii, 
suggesting a change of IR-bright galaxy fraction with 
radius. To quantify this trend, we plot the fraction of 
IR-bright galaxies at different projected radii in units of 
i?2oo- F° r this fraction calculation, we use the brighter R 
magnitude limit (Mr < —20.15). As shown in Fig. [71 the 
fractions in A3266 are similar to those in Coma, except 
for the point at R ~ 0.5i?2oo; where there are ~ 20% 
more IR-bright galaxies in A3266 than in Coma. 

Despite this difference, both clusters show a general 
increase of IR-bright galaxy fractions from the central to 
the outer region. To improve the statistics of this general 
trend, we combine the data from the two clusters and find 
that the fractions increase smoothly with radius, with 
a linear relation, (0.24 ± 0.02) + (0.36 ± 0.06)i?/i? 2 oo- 
This general trend is consistent with the observed de- 
crease of star-form ing galaxy fractions with density (e.g., 
iLewis et al.ll2002T ). In the Coma cluster, for which B 
band photometry is available, we also calculate the blue 
galaxy fraction as a function of radius. The blue galax- 
ies are defined as galaxies with B — R < 1.95 — 0.05.R, 
0.18 mag bluer than the mean color of the red sequence 
galaxies. Among cluster members with Mr < —20.15, 
the limit used to calculate the IR-bright galaxy fraction, 
the blue galaxy fraction increases from ~ 6% in the cen- 
tral region to ~ 15% at 0.7i?2oo- These fractions are 
4-5 times smaller than the IR-bright galaxy fractions. 
However, if we increase the lower IR luminosity cut to 
Lju > 10 43 - 35 ergs s^SFR > lM^r" 1 ), the IR-bright 
galaxy fractions and the blue galaxy fractions match each 
other quite well. The colors of those IR-bright galaxies 
with SFR > lAlQ-yr^ 1 confirm that the majority are 
blue galaxies. This indicates that the decrease of the 
blue galaxy fraction with galaxy density can only account 
for part of the decrease of the IR-bright galaxy fraction, 
namely, just for those galaxies with SFR > 1 M Q yr _1 , 
and that the fraction of the red galaxies with moderate 
SFR (up to 1M Q yr _1 ) also decreases with galaxy den- 
sity. 

4.5. The Evolution of the Integrated SFRs 

In iBai et all (|2007t ), we studied the integrated SFRs 
(ESFRs) of more than 10 clusters and found some tenta- 
tive correlations with redshifts and masses. To test those 
results, here we add the data from A3266. We integrate 
the SFRs of the cluster members within the 0.5i?2oo (~ 
1.5 Mpc) down to a SFR limit of 2 M yr' 1 . This gives a 
value of 39 M Q yr _1 . The spectroscopic incompleteness 
correction slightly increases this value to 43 M yr^ 1 . 
The radius (0.5i? 20 o) and the SFR limit (2 M yr^ 1 ) for 
calculating the ESFRs are chosen to accommodate most 
of the cluster data in the sample. In A3266, this inte- 
grated SFR only accounts for a quarter of the total inte- 
grated SFR of all the cluster galaxies that we study. We 
also calculate the mass-normalized ESFR (ESFR/M), 
for which we use a cluster mass of 3.3 ± 0.1 x 10 15 M Q . 



The cluster ma ss is deduced using a velocity dispersion 
of 1255 km s" 1 (IChristlein fc Zab ludofi 2003[) and utiliz- 
ing the formula given bv lFinn et alj (|2004h . The error on 
the mass takes into consideration the en hanced velocity 
dispe rsion caused by a possible merger ( Quintan a et al.l 
fl996h . 

As can be seen in Fig. U both ESFR and ESFR/M 
in A3266 are consistent with the evolution trend that we 
found before. On the other hand, the relation between 
the ESFR/M and cluster mass is less conclusive com- 
pared to the evolution trend. Especially for the less mas- 
sive clusters, there is large scatter in both the ESFR and 
the ESFR/M. The four rich clusters studied by MIPS, 
despite the large mass range (see Panel d of Fig. [5]) , show 
good agreement in ESFR/M at the same redshift. This 
argues against a strong dependence of ESFR/M on clus- 
ter mass for these rich clusters. Similarly, the IR LFs of 
these four rich clusters also argue against a strong de- 
pendence on cl uster mass. This is consistent with the 
recent results of iFinn et al.l (|2008l ). who suggest that the 
difference between the fraction of star forming galaxies 
in the field and clusters first appears in clusters with ve- 
locity dispersions lower than 600 km s _1 (ie, groups); 
clusters with higher velocity dispersions all have similar, 
and much lower than the field, fractions of star forming 
galaxies. 

Using only the data points of these four clusters, for 
which we are most confident, we derive the evolution in 
the trend ESFR/M oc (1 + z) 5 3±1 - 2 . This relatively 
well-constrained evolution trend deduced from the rich- 
est clusters at two different epochs, lies under a more 
scattered evolution traced by other data points from ISO 
observations and H a surveys. The scatter may be partly 
due to the uncertainties of the aperture correction, de- 
tection limit correction and extinction correcti on (for H a 
data) we made to obtain ESFR (see details in IBai et ah! 
[2007f l. but more likely is due to some intrinsic difference 
in cluster properties that may affect their star formation 
properties. For example, one of the two clusters deviat- 
ing most from the evolution trend, CL 1040 (z = 0.7) has 
a much smaller mass (~5x 10 13 M Q , iClowe et a"Lll2006f> 
compared to other clusters in the sample and therefore 
probably has properties more similar to galaxy groups 
than clusters. The other one, CL 0024 (z = 0.4), is X- 
ray underluminous compared to the clusters with similar 
mass, suggesting a lack of hot intracluster medium (ICM) 
that may relate to the abnormally high star formation 
activity. To disentangle all these effects from the evo- 
lution, we need not only a comparison of a well-defined 
and much larger cluster sample at different redshifts, but 
also comparison of cl uster and group s ample s. 

A recent study by ISaintonge et al.l (|2008f ) found that 
the fraction of star forming galaxies in eight massive clus- 
ters (M > 5 x 10 14 M Q ) increases from 3% at z = 0.02 to 
13% &tz = 0.83. The SFR limit cut (> 4 M Q yr" 1 ) and 
the region (<1 Mpc) they used to calculate the fractions 
are comparable to what we used to calculate the ESFR. 
The star-forming galaxy fraction evolution, ~ (1 + z) 2 ' 4 , 
is not as strong as the ESFR/M evolution we find here. 
This is because the evolution of ESFR/M also takes 
into account the evolution of L* IR for star-forming galax- 
ies. This again confirms what we find by studying the 
evolution of the IR LFs: both the density and the L/r 
of the star-forming galaxies have evolved strongly from 
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z = to z = 0.8 in rich clusters. 

5. DISCUSSION 

5.1. The Universality of the IR LF 

The overall IR LFs of A3266 and the Coma cluster 
are similar both in shape and in normalization at Ljr > 
10 42 ' 7 ergs s -1 . Since both of the observations cover 
about 70% of the i? 2 oo area of the cl usters, field coverage 
is not a factor in the c omparison (iPopesso et al.l 120061 : 
iBarkhouse et al.ll2007T ). iBai et all (|2006h compared the 
IR LF of the Coma cluster with that of the local field 
and did not find a significant difference in L* IR . The 
comparison of the faint end slope gives different results 
depending on which field IR LF is used, but generally no 
significant difference is found. Since the IR LF of A3266 
is similar to the Coma IR LF, it confirms the constancy 
of the shape of the IR LF from clusters to the field at 
the bright end (Lju > 10 43 ergs s _1 ), even for very rich 
clusters with M~ 10 15 M< 



ICortese et afl (|2003l I2005D studied the UV LF of sev- 
eral local clusters and found that the shape of the LF of 
the star-forming gal axies in clusters is co n sistent with 
that of the field. Ilglesias-Paramo et al.1 ((2002) con- 
structed the Ha LF of three local clusters and found 
it is consistent with the field Ha LF at the bright end. 
IChristlein fc Zabludofj (|2003h found that the R-band LF 
of star-forming galax ies in clusters is indist inguishable 
from that of the field. iDe Propris et alj (|2003[ ) also found 
that the B-band LF of the late-type galaxies in clusters 
is very similar to that of field galaxies. These results sug- 
gest a universal shape of the LF of star-forming galaxies 
in very different environments. 

A possible explanation for this universal shape is that 
the majority of the star-forming galaxies we see in the 
clusters are galaxies accreted from the field and their 
star formation has not been altered by the cluster envi- 
ronment yet. The accreted galaxies referred to here are 
not interloper galaxies whose redshifts happened to fall 
into the redshift range of the cluster but which are not 
physically bound to the cluster. The normalization of 
the cluster LF is still 10 times larger than the field IR 
LF normalization multiplied by the physical scale corre- 
sponding to the cluster redshift range, indicating the con- 
tamination from such interlopers is small. The fact that 
the velocities of the galaxies studied here fall within the 
maximum line-of-sight velocity requirement for a bound 
cluster member (see Fig. [IJ also confirms that they are 
not interlopers. 

Although these galaxies are not interlopers and their 
projected radial distances are all within -R20O1 it might 
be that they are recently accreted and still lie outside of 
R200, where some of the most important cluster environ- 
mental processes, e.g., ram-pressure stripping, tidal in- 
teraction, and starvation, have not started to affect them 
yet. If we assume the mass function of the cluster galax- 
ies does not change across the cluster and the cluster 
galaxy distribution follows the dark matt er distribution, 
which can be modeled by a NFW profile (jNavarro et al.l 
Il997f ) w ith i?2nn/'Ts = 6 for a cluster at the mass range of 
A3266 (|Eke et all 1 1991) . we find that about 30% of the 
galaxies projected within -R200 lie beyond i?2on- However , 
many simulations (jBalogh et all 120001 : iGill et al.l |2005| ) 
have shown that even for galaxies in the 1-2 -R200 ra- 
dius, about 50% are a 'backsplash' population, galaxies 



which have once penetrated into the -R200 region. This 
implies that only about 15% of the cluster members 
we studied here could really be recently accreted field 
galaxies that have never been inside of i?2oo- Therefore, 
even if all of the recently accreted field galaxies are star- 
forming galaxies, they cannot entirely account for the 
35% of star-forming galaxies we detect in A3266 down to 
M R < -19.15 (23% in Coma). 

Given this evidence, we conclude that hot gas starva- 
tion is unlikely to alter the overall star formation proper- 
ties of cluster galaxies. This is because starvation starts 
to work around i?200 and it is a slow process (~Gyr). 
If it does have a major effect on galaxy star formation, 
we would expect to see many star-forming galaxies that 
are in the transition stage of the process and they should 
have a lower SFR (and L* IR ) compared to the field galax- 
ies not affected by this process. This contradicts the un- 
changed bright end IR LF shape we find from the field 
to the cluster. However, we need to note that given the 
accuracy of our measurement of the cluster IR LF, we 
are not sensitive to a variation of 30% in L* IR . In addi- 
tion, the uncertainty of the field IR LFs further decrease 
the sensitivity of the comparison of environments. So 
a modest change of L* IR value between cluster and field 
galaxies might be buried by the uncertainties. 

Nevertheless, our result is compatible with certain 
classes of cluster environmental effects. For example, if 
an effect is limited to a smaller central region than i?200i 
it would be possible that these star forming galaxies are 
cluster members that have not penetrated deep inside 
the cluster, and therefore their star formation remains 
unaffected. For instance, ram-pressure stripping only 
works efficie ntly within 0.5 -R 200 for a Milky Way- type 
galaxy (e.g.. lTreu et al.ll2003l ). typical in stellar mass of 
the galaxies of the star-forming galaxies in our study. Us- 
ing the same NFW model, we can estimate that about 
60% of the galaxies with projected radial distance < -R200 
lie beyond 0.5i?2O0- Of the backsplash population in the 
1-2 i?2oo radiu s, about 90% h ave penetrated into the 
0.5i?2oo region (IGill et all I2005h . If we assume galaxies 
located between 0.5i?2oo and -R200 have the same frac- 
tion of the deep-penetrating backsplash population, we 
find about 30% of galaxies with projected radial distance 
< -R200 have never been inside of the 0.5i?2oo region. 
This fraction is comparable to the star-forming galaxy 
fraction we find. If we take into consideration that the 
deep-penetrating fraction in the 0.5-1 -R200 radius range 
may be higher than in the 1-2 i? 2 oo radius range and 
the fact that not all field galaxies are star-forming galax- 
ies to start with, then we still fall a little short to fully 
account for the star-forming galaxies detected in the clus- 
ter. However, because the effective ram-pressure strip- 
ping radius c an be less than .5i?2oo for galaxies not in 
radial orbits (jTreu et al.ll2003l ). the percentage of galax- 
ies not affected by ram-pressure stripping can be larger 
than 30% and therefore we can account for the majority 
of the star-forming galaxies we find in the cluster. 

In addition, because the ram-press ure stripping works 
on a short time scale (~ 10 7 yrs, lAbadi et alj 119991 ) 
and can convert a star-forming galaxy into a quiescent 
one very quickly, th e re will be few galaxi es in transition 
(|Balogh et al.l l2004t IVerdugo et al.l l2008h and this will 
keep the shape of the IR LF unchanged while changing 
the fraction of star- forming galaxies. To keep the shape 
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of IR LF unchanged, the ram-pressure stripping also 
needs to work equall y effectively on the relevant range of 
galaxy luminosities (jCortese et al.ll2003h . This seems in- 
consistent with our understanding of ram-pressure strip- 
ping, which should be more efficient for the low-mass 
galaxies than for the high-mass ones. But in the Coma 
cluster where we can constrain the faint end slope of the 
IR LF, we do see a flatter faint end slope in the core re- 
gion, which may be due to the more efficient ram-pressure 
stripping for fainter galaxies. 

As a conclusion, our results suggest that the cluster 
environment has little effect on the shape of the high- 
end cluster infrared luminosity function, but is charac- 
terized by a decrease of IR-bright galaxy fraction towards 
the cluster center. Nonetheless, the proportion of such 
galaxies is high enough to make it implausible that they 
are all recently accreted from the field. To be compatible 
with both the number of infrared-active galaxies and the 
behavior of the LF, any significant cluster environmental 
effect must be largely confined to the inner region (e.g., 
0.5 -R200) an d probably acts relatively quickly. These ar- 
guments would be compatible with effects through ram- 
pressure stripping or a similar process, but not with star- 
vation. 

5.2. The Universality of the Evolution of the IR LFs 

In addition to the universality of the shape of the 
bright-end IR LF from clusters to field at z = 0, the 
evolution of this LF also shows universality from clus- 
ters to field. If this evolution is passive or driven by 
some mechanism that works in both environments, then 
it is no surprise that we see similar evolution in clusters 
and field. However, our observations are also compatible 
with some specific forms of differential evolution in the 
cluster environment that might arise from cluster-specific 
mechanisms such as ram pressure stripping. In general, 
slow mechanisms are unlikely because they would leave a 
transitional galaxy population in the process of evolving, 
and we would have to appeal to coincidence to explain 
the similarity in L* IR between the field and the clusters. 
However, as discussed in the preceding section, a fast- 
acting cluster environmental effect is possible, if the clus- 
ter keeps replenishing its star forming galaxies from the 
field. In this case the cluster mechanism can quickly turn 
off the star formation in some of accreted galaxies, for 
example, those penetrate deep inside the cluster, leaving 
the L* IR of rest of the star forming galaxies untouched. 
As a result, the observed L* IR evolution in cluster will 
just follow the Lj R evolution in the field, which, on the 
other hand, can be driven by some environmental mech- 
anism not efficient in the cluster, e.g., galaxy merger. 

The interpretation of the evolution of the IR LF nor- 
malization, however, is more complicated. This normal- 
ization depends both on the accretion history and SF- 
suppression efficiency of the cluster. It also depends on 
the space density of the star-forming galaxies in the field, 
which is the reservoir for cluster accretion. For example, 
assume rich clusters at high redshift have higher accre- 
tion rates than local clusters, as suggested by some hi- 
erarc hical structure formation studies (e.g.g iKauffmannl 
1995), and also have comparable SF-suppression effi- 
ciency to local clusters. Then the evolution of the IR LF 
normalization would be stronger than that in the field, 
but would still more or less follow the field number den- 



sity evolution. This behavior would be compatible with 
our results. However, given the large uncertainties in 
both the normalization evolution and the relevant pa- 
rameters of high-z cluster behavior, this interpretation is 
very uncertain. 

6. CONCLUSIONS 

Despite its complicated dynamical state, the IR LF 
in A3266 is very similar to that of the Coma Cluster, 
suggesting a universal form of the bright end IR LF 
(Lin > 10 43 ergs s _1 ) in nearby rich clusters. The 
shape of this cluster IR LF, at least at the bright end, 
is not significantly different from the field IR LF at the 
same redshift. The fraction of the IR-bright galaxies with 
Lm > lO 42 ' 7 ergs s' 1 (SFR> 0.2 M Q yr~ l ) down to 
Mr < -20.15 is 40±|% in A3266 and 31+5% in Coma. 
There are few IR-bright galaxies within a projected ra- 
dius of r < 0.3 Mpc, and even they may not actually 
lie within the cores. The fraction of IR-bright galaxies 
in these two clusters increases linearly with radius, as 
(0.24 ± 0.02) + (0.36 ± 0.06)i?/i? 2 oo- The decrease of 
the blue galaxy fraction toward higher density regions 
only accounts for part of the trend, specifically, for the 
IR-bright galaxies with SFR> IMq yr _1 ; the fraction of 
red galaxies with moderate SFR (0.2 M Q yr^ 1 < SFR < 
1 Mq yr^ 1 ) also decreases with increasing galaxy den- 
sity. These results show that the cluster environment is 
characterized by a decrease of IR-bright galaxy fraction 
toward the cluster center, but not by a change in L* IR . 
If some cluster environmental mechanism is responsible 
for the difference, this mechanism must be confined to 
the central regions (e.g., < 0.5i?2oo) of the cluster and 
probably acts relatively quickly. These properties would 
be compatible with a process like ram pressure stripping, 
but not with a slower-acting one such as starvation. 

When comparing the A3266-Coma composite IR LF 
with that of two rich clusters at z ~ 0.8, MS 1054 and 
RX JO 152, we find a strong evolution characterized by 

L* IR oc (1 + zf- 2 -", $* m oc (1 + z) 1 - 7 - 1 !-". This L* IR evo- 
lution is indistinguishable from that in the field, and the 
evolution is stronger, but still consistent with that in 
the field. The similarity of the evolution of bright-end IR 
LF in very different cluster and field environments sug- 
gests either this evolution is driven by the mechanism 
that works in both environments, or clusters continually 
replenish their star-forming galaxies from the field, yield- 
ing an evolution in the IR LF that is similar to the field. 

The mass-normalized integrated SFR (ESFR/M) 
within 0.5i?2oo (~ 1-5 Mpc) of A3266 and Coma are sim- 
ilar and show an evolution of (1 + z ) 5 - 3±1 - 2 compared 
with those of the two high-z clusters. The evolution de- 
rived from these four rich clusters lies under a more scat- 
tered evolution trend traced by a larger cluster sample. 
The large scatter indicates that in addition to evolution, 
cluster star formation may correlate with cluster masses, 
IGM properties, and perhaps other parameters. How- 
ever, for the four rich clusters studied by MIPS, neither 
the IR LFs nor "ESFR/M within 0.5i?2oo show a depen- 
dence on cluster mass, suggesting the mass dependence, 
if any, is weak for clusters with M >~ 4 x 10 14 Mq. 
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Fig. 1. — The line-of-sight velocity relative to the cluster mean versus projected distance from cluster center. The open circles are cluster 
members selected by ±3x velocity dispersion, and the filled ones are those detected at 24 /im. The symbol sizes are proportional to the 24 
(im flu x. The two dashed curves are the most stringent criterion of the maximum line-of-sight velocity proposed by den Hartog & KatgerU 
(1996). All the members fall inside (or very close) to these curves. The three vertical lines at 0.13, 0.3 and 0.55 i?200 are the boundaries 
of the three regions used in §4.3 and §4.4. 
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Fig. 2. — The completeness of the spectroscopic data in the cluster field as a function of R magnitude. The error bars are the lc errors 
given by binomial statistics. 
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Fig. 3. — The IF LF of A3266 (open squares). The filled squares are the IR LF after correction for the incompleteness of t he spectroscopic 
data and the solid curve is its best fitting Schechter function. The gray filled circles are the IR LF of the Coma cluster I Bai et al]|2006l ) 
and the dash-dotted curve is its best fitting Schechter function. The dotted vertical line is the luminosity corresponding to the detection 
limit of the 24 /im data. The dashed vertical line is the luminosity corresponding to the spectroscopic detection limit. 




Fig. 4. — The composite IR LF of A3266 and Coma (filled stars). The best fitted Schechter function is shown as the solid curve. 
The filled circles are the composite IR LFs of MS 1054-03 (z = 0.83) and RX J0152 (z = 0.84). The dash-dotted line is its best fitted 

Schechter function. The dotted line is the IR LF of A3266 and Coma evolved to z = 0.83 using L* R oc (1 + zf -0.2 ; « (1 + z) v -a.s 

the Floc'h et alj|2005h . 
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Fig. 5. — The sky map of A3266 and Coma. East is to the left and north is at the top. The IR cluster members are shown as filled 
circles. The symbol sizes are proportional to the 24 (im flux. The open circles are spectroscopically confirmed cluster members. The three 
big circles indicate the three regions within which we extracted IR LFs (r = 0.13, 0.3, 0.55il20o). 
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Fig. 6. — The composite IR LFs in different regions of the clusters. The solid curves are the best fitted Schechter function of the total 
composite IR LF of these two clusters. The dotted curves are the best fitted Schechter function in the different regions. 
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Fig. 7. — The fraction of star-forming galaxies (logLm > 42.7, ~ SFR > 0.2 Mq yr~ l ) as a function of projected distance from the 
cluster centers. The fractions are only calculated for galaxies with Mr < —20.15. The error bars are the la errors given by binomial 
statistics. 
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Fig. 8. — (a) and (6), the integrated SFRs (< 0.5iJ20o) vs. rcdshifts and cluster masses; (c) and (d), the mass-normalized SFRs 
(< 0.5i?20o) vs - redshifts and cluster masses. The filled stars are four clusters observed with MIPS: Coma, A3266, MS 1054-03, and 
RX J0152. Filled circl e s are the clusters observed with ISOCAM, and open circles are from Ha emission line measurements, adjusted as 
described in Bai et al. (2007). The dotted curve in (c) is the fitted correlation between mass-normalized SFRs and redshifts for the four 
clusters observed by MIPS. 
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Note. — Table[T]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 

regarding its form and content. 

References for redshifts: ( 1) IChristlein fc Zabludofll (2003); (2) IQuintana et al.l (Tl996| s ); QI IGreen eTaLl fl990). 



